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We calculate the cross section of a pion diffraction dissociation in two jets with large trans- 
verse momenta originating from a hard gluon exchange between the pion constituents. To 
the leading logarithmic accuracy (in energy), the contribution coming from small trans- 
verse separations between the quark and the antiquark in the pion acquires the expected 
factorized form, the longitudinal momentum distribution of the jets being proportional 
to the pion distribution amplitude. The hard gluon exchange can in this case be consid- 
ered as a part of the unintegrated gluon distribution. Beyond the leading logarithms (in 
energy) this proportionality does not hold. Moreover, the collinear factorization appears 
to be broken by the end-point singularities. Remarkably enough, the longitudinal mo- 
mentum distribution of the jets for the non-factorizable contribution is calculable, and 
turns out to be the same as for the factorizable contribution with the asymptotic pion 
distribution amplitude. 



1. It has been conjectured a long time ago |1|] that pion diffraction dissociation 
on a heavy nucleus tcA XA is sensitive to small transverse size configurations of pion 
constituents. It was later argued H, Q that selecting a specific hadronic final state that 
consists of a pair of (quark) jets with large transverse momentum gi^ ~ —q2i. one can 
obtain important insight into the pion structure as it turns out that the longitudinal 
momentum fraction distribution of the jets follows that of the pion valence parton con- 
stituents. A measurement of hard dijet coherent production on nuclei presents, therefore, 
the exciting possibility of a direct measurement of the pion distribution amplitude and 
provides striking evidence @] that this distribution is close to its asymptotic form. 

From the theoretical point of view, the principal question is whether the relevant 
transverse size of the pion r±_ (alias the scale of the pion distribution amplitude /i = 1/t_l) 
is determined by the color transparency condition /i ~ A^/^Aqcd or whether it is of the 
order of the transverse momenta of the jets fi ^ q±. In the latter case one could envisage 
a factorization formula for the amplitude of hard dijet coherent production of the type 
(cf. Fig. P 
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M^^2iets = J dz j dxi(f)^{z',nl)TH{z',Xi,n%) F^{xi,nl) . (1) 



Here (f).„{z' , fip) is the pion distribution amplitude, and F^{xi,^p) is the non-forward 
(skewed) gluon distribution P, || in the target nucleon or nucleus; Xi and X2 = Xi — ( 
are the momentum fractions of the emitted and the absorbed gluons, respectively. (The 
asymmetry parameter ( is fixed by the process kinematics, see below.) Th{z' ,xi, fi'jp) is 
the hard scattering amplitude and /ii? is the (collinear) factorization scale. 

To the best of our knowledge, this question has never been studied. As a first step 
in this direction, in this letter we present an explicit calculation of the leading-order 
contribution to the imaginary part of Th{z', Xi, n'jp) corresponding to a single hard gluon 
exchange. Since we are not interested in the color transparency phenomena, we restrict 
ourselves to scattering from a single nucleon. We will find that the structure of the hard 
gluon exchange is such that it generates an enhancement by a logarithm of the energy 
in the region ~ z where z is the energy fraction carried by the (quark) jet. If only 
this logarithmic contribution is retained, the longitudinal momentum distribution of the 
jets indeed follows the shape of the pion distribution amplitude [@, 0. The hard gluon 
exchange can in this case be considered as a part of the unintegrated gluon distribution, 
as advocated in [Q. Beyond the leading logarithms (in energy) this proportionality does 
no longer hold. Moreover, the collinear factorization (^ appears to be broken by the 
end-point singularities. Remarkably enough, the longitudinal momentum distribution of 
the jets for the non-factorizable contribution is calculable, and turns out to be the same 
as for the factorizable contribution with the asymptotic pion distribution amplitude. 

2. The kinematics of the process is shown in Fig. |l]. For definiteness, we consider 
pion scattering from a nucleon target. The momenta of the incoming pion, incoming 
nucleon and the outgoing nucleon are pi,P2 and respectively. The pion and the 
nucleon masses are both neglected, pf = 0, pi = {p'2)'^ = 0. We denote the momenta of 
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Figure 1: Kinematics of the colierent hard dijet production tt 2 jets. The hard scattering 
amplitude Th contains at least one hard gluon exchange. 

the outgoing quark and antiquark (jets) as qi and q2, respectively. They are on the mass 
shell, qf = = 0. We will use the Sudakov decomposition of 4-vectors with respect 
to the momenta of the incoming particles pi and p2- For instance, the jet momenta are 
decomposed according to 

qi = zpi + —P2 + gi± , q2 = zpi + —P2 + g2± (2) 

zs zs 

such that z is the longitudinal momentum fraction of the quark jet in the lab frame. 
We will often use the shorthand notation: u = (1 — u) for any longitudinal momentum 
fraction u. The Dirac spinors for the quark and the antiquark are denoted by u{qi) and 
v{q2)- 

We are interested in the forward limit, when the transferred momentum t = {p2 —p'2)^ 
is equal to zeroQ, and the transverse momenta of jets compensate each other qi± = q±, 
q2± = —q±- In this kinematics the invariant mass of the produced qq pair is equal to 
= q\/zz. The invariant cm. energy s = {pi + ^2)^ = '^PiP2 is taken to be much 
larger than the transverse jet momentum q±. In what follows we neglect contributions 
to the amplitude that are suppressed by powers of 1/s. 

The general scheme of the calculation can be explained as follows. Since the hard 
scattering amplitude Th{z',xi, fij;,), by assumption, does not depend on the target, we 
choose to consider the hard dijet production from a quark, nq {qq)q- We replace the 
pion by a collinear qq pair with the momenta z'pi and z'pi, respectively. The probability 
amplitude to find a particular value of the momentum fraction z' is given by the pion 
distribution amplitude defined as 

1 

{0\d{yh,^,u{-y)\7c+{p))y.^o = ip, U j dz' e'^^^'-'^^^y^ (t^^{z') , (3) 



*If the target mass m is taken into account, the momentum transfer t = [p2 — contains 
a non- vanishing longitudinal contribution and is constrained from below by \t\ > tQ, where to = 
{inn?M'^)/[s — TTn?)"^ ^ A'P being the invariant mass of the dijet. 
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(a) 



Figure 2: Cut diagrams (examples) for the imaginary part of the amplitude vrg {Qq)Q- The 
cut quark (gluon) propagators are indicated by crosses. 



where /^r ~ 131 MeV is the pion decay constant. Light-cone dominance for the t-channel 
gluon emission is not assumed from the beginning, but has to follow from the calculation 
in order that the result can be interpreted in the sense of the factorization formula 
(HI). To this end, the gluon transverse momentum k± is kept nonzero and we show 
that the amplitude in question contains a collinear logarithm lng^//i|. coming from the 
integration region fip kj_ <^ q"^. This property allows to calculate the upper part of 
the diagrams represented schematically in Fig. |I| assuming that k'j_ <C i.e. in the light- 
cone limit. In the last step, the collinear logarithm is interpreted as a contribution to the 
gluon distribution in the target and the (perturbative) non-forward gluon distribution of 
the quark is substituted by the (nonperturbative) non-forward gluon distribution in the 
nucleon. 

In this letter we calculate the imaginary part of the amplitude. The real part can in 
principle be restored from the imaginary part using dispersion relations. One rationale 
for this procedure is that at high energies the scattering amplitudes corresponding to 
Pomeron exchange are dominated by their imaginary parts. The second rationale is 
simplicity: the corresponding cut diagrams (see Fig. ^ for examples) are built of tree-level 
on-shell scattering amplitudes and their form is strongly constrained by gauge invariance, 
see below. 

The existing cut diagrams can be grouped into the four gauge-invariant contributions 
shown in Fig. ^-d, which differ by the position of the hard gluon that provides the large 




(a) (6) (c) id) 

Figure 3: The decomposition of the imaginary part of the amplitude nq {QQ)q into four 
gauge-invariant contributions. 



momentum transfer to the jets. The corresponding contributions to the amplitude will 
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be denoted as M{a), -^(c) and J^[d)- For example, in Fig. |^a it is assumed that the 

hard gluon exchange appears to the left of the cut. This contribution is given by the sum 
of 10 Feynman diagrams one of which is shown in Fig. Similarly, the contribution in 
Fig. is given by the sum of 10 diagrams with the hard gluon exchange appearing to the 
right of the cut; a typical diagram is shown in Fig. The two remaining contributions 
in Fig. 1^ and Fig. ^ take into account the possibility of real gluon emission in the 
intermediate state. The filled circles stand for the effective vertices describing the gluon 
radiation, see Fig. ^. Each of the two contributions in Fig. ^,d corresponds to a sum of 
9 different Feynman diagrams, see Fig. ^ for an example. 
We use Feynman gauge and perform the usual substitution 

in the propagators of t— channel gluons. This is correct up to terms ~ 1/s which is 
sufficient for our purposes. The t— channel gluon momenta are expanded according to 

k'^ = aip^ + xi P2 + k^, 

= + a;2P2 + ^± • (5) 

As is easily checked by inspection, in any cut diagram two internal lines are on the 
mass shell. The corresponding two on-shellness conditions fix ai and xi and relate the 
variables z' and xi to one other. In addition, 02 and X2 are fixed by the energy-momentum 
conservation. Since a(i,2) and a;(i^2) are all of the order of 1/s, the 1/A;^2 factors in the 
propagators of the t— channel gluons can be approximated by k"^^ 2) = tt(i,2)a^(i,2)S — ^1 = 
—k"^ + 0{l/s). Using the on-shellness conditions for the contributions in Fig. |^ and 
Fig. ||b one obtains Xi = ^2 = 0, for any z'. For Fig. |^d one finds Xi = (z'z/{z' — z), 
X2 = C,zz' I {z' — z) and z' > z, where the last condition ensures that the energy of the cut 
gluon is positive. Finally, for the set of cut-diagrams corresponding to Fig. |^ we obtain 
Xi = (zz'/{z — z'), X2 = Cz'z/{z — z') and z > z' . 

After the on-shellness conditions are used, a single integration over the gluons trans- 
verse momentum k^ remains: 

Jup{k±, q±) Jdown{k±, q±) , (6) 

where fc^ comes from the product of the two gluon propagators. J„p and Jdown are dubbed 
impact factors and stand for the upper and the lower parts of the diagrams Fig. |a-d 
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(connected by the two-gluon exchange). The representation (|^) is well known |^ from 
QED scattering at high energies. 

Properties of the impact-factors Jup and Jdown as a functions of k±_ at fc_L — > are 
of crucial importance. Since Jdown is the impact-factor of a point-like target quark, 
Jdown{k±,q±) ~ const. On the other hand, Jup{k±,q±) stands for the scattering of the 
colorless qq (Fig. ^-b) or qqG (Fig. 0c-d) state having a transverse size ~ l/g_L and has 
to vanish at small k± -C g_L, Jup{kj_, q±) ~ k\, as a consequence of the color neutrahty of 
the quark-antiquark pair: A gluon with a large wave length ~ l/fc^ cannot resolve a color 
dipole of the small size ~ l/q±- Since in our case there are two gluons, Jup is proportional 
to the product k±-k± = k'jj\. In the opposite hmit of large transferred momenta, A;^ ^ q±, 
the two t— channel gluons are forced to couple to the same parton (quark or gluon) in 
the upper block in Fig. ^-d. It follows that at large k± Jup{k±, q±) ~ const. 

Taking into account the above properties of the impact-factors we conclude that 
the transverse momentum integration in (^ diverges logarithmically at small k_i_ and 

2 

the integral can be estimated as ~ J"-^ dk'j_/k'j_ ~ Ing^., as expected. The region of 
k'j_ > q\ does not produce the large logarithm and can be neglected. Note that the correct 
small k± behavior of the impact factors is only recovered in the sum of cut diagrams for 
the gauge invariant amplitudes M.{a)i -^(fe), -^(c) and M.{d)i but not for each diagram 
separately. 

In addition to the diagrams discussed so far, the amplitude nq {qq)q receives a 
contribution from the three-gluon exchange in the t-channel. Such terms can be viewed 
as belonging to the cut diagrams shown in Fig. ^ in which the hard gluon in the blob 
is attached to the bottom quark line. We have checked that this extra contribution does 
not contain the large coUinear logarithm Ing^. and therefore we neglect it. 

3. To start with, consider the calculation of J^[d)- Let /'^ = aiPi + xiP2 + be the 
momentum of the (real) gluon in the intermediate state and let e^{l) be one of the two 
physical polarization vectors. The two conditions (e ■P2) = and (e ■ /) = fix the gauge 
and result in e^{l) = e'^ + {e±l±) /{ai s). 

The effective vertex corresponding to the sum of the three diagrams in Fig. |^ has the 
form 



^ 1 I 



1 (t4")..-- itt').. 



u{qi 



s 



u{z'pi) . (7) 



Here and are the SU (3) generators. The color indices / and a belong to the emitted 
gluon and the t— channel gluon, respectively. We have also introduced an auxiliary two- 
dimensional vector 6^ defined as: 

= _2i^^g^, (8) 

Note that the effective vertex, in the limit of small k±, is proportional to b cc k±. The 
constant terms cancel in the gauge invariant sum of the diagrams in Fig. ^. 

^The 0{k\) behavior can be traced to the gauge invariance of the amphtude, see @] for the details. 
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The second effective vertex in Fig. ^ has a similar form. Combining both of them 
and performing the sum over the polarizations of the emitted gluon we obtain the impact- 



factor Jnp . Since each effective vertex is proportional to k±, it follows that Jup' ~ k'^, as 
expected. The result for the amplitude M{d) is obtained using the representation in (H). 
The calculation of A^(c) is very similar. The result for their sum reads: 



where 



id) 



DCl 



dz' (pTr{z') 



z z 



+ 1 X 



X 



z z 



+ 1 + 



1 



{N^ - 1) V^' z' 



+ 



Q(z' -z) ^ Q{z-z') 



[z' - z) 



[z - z' 



D = -isUa^ u{qi)'y5—v{q2) Oij dec' , 

ql s 



(9) 



(10) 



and Cf = {N^ — i)/'^Nc. The color indices correspond to the produced quark- 

antiquark pair (jets) and (c, c') stand for the color indices of the target quark in the 
initial and the final state. The contributions ~ Q{z' — z) and ~ 0(z — z') belong to Ai{d) 
and A^(c), respectively. 

For the cut diagrams in Fig. ^a and Fig. §3 we present the final results: 



M 



(a) 



M 



-DCl 



dk] 
dkl 



dz' 



z z 
zz{ — + — 



+ 



{m - 1) \z'z' 



zz 



+ 1 



'ir 



The transverse momentum integrals J dk^/k^ ~ Ing^ in (|]) and (|TT|) can be identified 
with the (perturbative) non-forward gluon distributions of a quark: 



— Up 

IT 



dk] 



a 



— Cf Ingl 



TT 



(12) 



To justify this substitution, note that qF^{x) can be calculated to first order in pertur- 
bation theory from the evolution equation p| 



9± 



2^ 



PI\x,z)=Cf 



I dzP^%x,z),F^{z,ql) 

X 

X ~ ( 



+ 1 



(13) 



Using qF^{z) = 6{1 — z) and taking into account that in the high-energy region (,x 
the quark-gluon kernel simplifies to P^'^{x, 1) = 2Cf, we arrive at the substitution rule in 

Hn our calculation Fi^{x) enters in the DGLAP region a; > C only, see M for more details. 
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Eq. (0). The final step is to replace qF^{x) by the (nonperturbative) non- forward gluon 
distribution in the nucleon which is a nontrivial function of the parameters [S: 



{N{p')\y,y^GUO)GUy)\N{p))y2^o = u{p') Mp) 



(yp) 



dx [e 



-ix{py) 



(14) 

Here G^^, is the gluon field strength tensor. u{p) and u{p') are the spinors of the initial 
and final nucleons. 

Our final result for the imaginary part of the amplitude for dijet production from a 
nucleon reads 

lmA1 = — zs/^a^ 



/v2 .^(gl)75-^(g2)X5^J 



(15) 



with 



J 



dz' 4>n{z', jJ^) 



zz 



z z 

- + -) + 

z' z' ■ 



2Nr z'z' 



zz 



+ 1 



+ 



zz 



+ 1 



zz 



+ 1 + 



2N^ \z' 



z z 

^ + ^ 



e{z' - z) 



z' -z 



(z'z 
z' — z 



Q{z 



[z - z' 



Cz'z 



(16) 



The differential cross section summed over the polarizations and the color of quark jets 
is given by 



ixi- 



dqidtdz 8N^ql'^' ' ^^^^ 

The factorization scale /x^ has to be of order of the transverse momentum of the exchanged 
gluon. The expression in Eq. (p!6|) presents the main result of this paper. 



4. The integrand in ( |16|) is singular at z' = z, i.e. when the longitudinal momentum 
fraction carried by the quark coincides with that of the quark jet in the final state, 
and at the end-points of the integration region z' — > and z' —>■ 1. Let us discuss the 
contributions from these regions in some detail. 

The singularity at z' = z is present in the contributions in Fig. ^c,d which include real 
gluon emission in the intermediate state. The logarithmic integral J dz' /\z — z'\ ~ Ins is 
nothing but the usual energy logarithm that accompanies each extra gluon in the gluon 
ladder. Its appearance is due to the fact the the gluon in Fig. ^,d can be emitted in 
a broad rapidity interval and is not constrained to the pion fragmentation region. To 
logarithmic accuracy we can simplify the integrand in (pIB]) by assuming z' = z everywhere 
except for the diverging denominators and the argument of the gluon distribution, to get 



X 



AN,<p^{z) 



dz' 



z' z 



z — z 



m- — Ai)^^N,<PM —Fdy.qi)- (is) 



dy 
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For a flat gluon distribution -F^(y) ~ const at ?/ — 0, and the integration gives const-ln 1/C, 
which is the above mentioned logarithm. Note that the color factors combine to produce 
Ca = Nc signaling that the relevant Feynman diagrams in Fig. ^,d are those with a 
three-gluon coupling. Moreover, the factor 2Nc/y appearing in ([TsD can be interpreted 
as the relevant limit of the DGLAP splitting function 

1 1 
^l^i^c(^ = C,gl) = ^/ dyP^^{(,y)F,{y,ql)c^l^ J dy^F,{y,ql). (19) 

The quantity on the l.h.s. of ([T9|) defines what can be called the unintegrated non-forward 
gluon distribution and the physical meaning of Eqs. (|18D and ([l9|) is that in the region 
z' z hard gluon exchange can be viewed as a large transverse momentum part of the 
gluon distribution in the proton, cf. This contribution is proportional to the pion 
distribution amplitude (j)T^{z,qj_) and contains the enhancement factor lnl/( ~ Ins/g^. 

Next, consider the contribution to the imaginary part of the amplitude of the dijet 
production coming from the end-points z' and z' ^ 1. Summing both of them and 
using the symmetry of the pion distribution amplitude we obtain 



Since 4>Tr{z') ~ z' at z' 0, the integral over z' diverges logarithmically. This diver- 
gence indicates that the coUinear factorization conjectured in (|lD is generally not valid. 
Remarkably, the divergent integral containing the pion distribution amplitude is just a 
constant and does not involve any z-dependence. Therefore, the longitudinal momen- 
tum distribution of the jets in the nonfactorizable contribution is calculable and, as it 
turns out, has the shape of the asymptotic pion distribution amplitude 4''^^{z) = Qzz. 
The corresponding physical process is the following. The limit z' ^ 1 corresponds to 
a kinematics in which the quark carries the entire momentum of the pion. The fast 
quark radiates a hard gluon which carries the fraction {1 — z) of quark momentum. 
This radiation is perturbative and is described by the effective vertex (|^) at 2;' = 1. 
At the final step the hard gluon transfers its entire longitudinal and transverse mo- 
mentum to the quark jet, and emits a soft antiquark which interacts nonperturbatively 
with the target proton and the pion remnant. The technical reason for the singular- 
ity at z' —>■ 1 can be traced to the quark propagator whose denominator has the form 
~ l/[qlz^/{z{z - z')) - 2z'{k_Lq_L)/{z - z') + k\{\ + z' jiz - z')) + ie], where is the 
transverse momentum of the t-channel gluon. In order to extract the leading small-Zc^ 
behavior corresponding to the logarithmic coUinear divergence we expand this denomi- 
nator at small k^. The leading term in this expansion cancels in the gauge invariant sum 
of diagrams, and the second term (linear in k^ ) plus similar terms from the nominator 
produce 1/(1 — z'Y' ■ It follows that the divergent logarithm in J /^^^{z') / z'"^ is of the form 
lng^//Xjj^ where /ijR is related to the average transverse momentum of the quarks inside 
the pion. It is possible that in the case of scattering from a heavy nucleus /xir may grow 
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as ~ because of color filtering. A detailed discussion of this effect goes beyond the 
tasks of this letter. 



5. In order to make an estimate we assume that skewedness does not affect signifi- 
cantly the X dependence and adopt the following simple model for the non-forward gluon 
distribution: 



Fc(x,/i) = F^(C,/i) 



(21) 



with A = 0.3 ||TT[. We further notice that at high energies, or at small there exists 
a mechanism which effectively suppresses the nonperturbative end-point contributions, 
or, better to say, enhances the factorizable contributions in comparison to the end-point 
ones. It is known that in the region of small x the scale dependence of the gluon density 
g{x,fi'^) is quite large and can be parametrized by the effective exponent 7: g{x,iJ?) = 
g{x, ^Di^fi"^ / fiiy . We take 7 = 0.3 and assume the same scale dependence for the 
non-forward distribution: 



FdC,^^') = FdC^^t>)(-2 ■ (22) 

It follows that the contributions coming from larger scales (alias smaller transverse dis- 
tances) are enhanced in comparison with the contributions coming from smaller scales 
(larger transverse distances). The typical virtualities of the quark and gluon propagators 
in the hard subprocess are ~ q'j_ when z' ~ z. They decrease when z' is close to the end 
points: ~ {z' / z)q'j_ at 2;' — >• and ~ {z'/z)q'j_ at z' 1, as the reminder that the soft 
contribution originates from large transverse distances. Therefore, it is natural to choose 
as the factorization scale: 

f.' = ql — . (23) 

zz 

With this choice we observe that the end-point singularities (formally) disappear. 

The longitudinal momentum fraction dependence of the jets calculated using the 
expression in ( |T6l) and the model for the nonforward gluon distribution in (|2T|) , (^2]) is 
shown in Fig. ^for two different choices of the pion distribution amplitude: (p'^iz) = 6zz 
and ^^^(2) = 30zz{2z - if P[. For this plot we have taken s = 1000 GeV^ and 



q± = 2 GeV which roughly corresponds to the kinematics of the E791 experiment 0]. 
The overall normalization is arbitrary. Notice that the Chernyak-Zhitnitsky (CZ) ansatz 
for the pion distribution amplitude leads to a much larger (integrated) cross section than 
the asymptotic distribution]^ This signals that the leading logarithmic approximation 



in Eq. (|T8D is not sufficient and regions of integration in (|T6D other than z ^ z' play 
an important role. One sees that the sensitivity to the shape of the pion distribution 
amplitude remains, although there is no direct proportionality as it was assumed in 

To summarize, in this letter we have calculated the hard gluon exchange contribution 
to the imaginary part of the amplitude of pion diffraction dissociation to two jets with 



§The calculation shown in Fig. ^ serves the illustrative purposes only. In the data analysis one has to 
take into account that the CZ model [|lOj is formulated at a low scale and has to be evolved to /Lt ~ g_L 
for a meaningful comparison. 
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Figure 5: The longitudinal momentum fraction dependence of the jets for two different choices 
of the pion distribution amplitude: i'^i^z) = ^zz and = 30zz{2z — 1)"^ shown by the 

solid and the dashed curve, respectively. 

large transverse momenta. The answer is given in Eqs. (|l6l), (IT^) and its main feature is 
that the jet longitudinal momentum distribution is not simply proportional to the pion 
distribution amplitude. The actual dependence is rather elaborate and it has to be taken 
into account in the data analysis. Our result can be improved by calculating the quark 
contribution which may be important in the energy range of the E791 experiment and 
by elaborating on possible nuclear effects that were not taken into account in the present 
study. 
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Note added 



Our result in Eq. (|T^) does not agree with an independent calculation using a different 
method where the light-cone dominance was assumed from the beginning. We plan to 
investigate the reasons for this disagreement in a separate publication. 
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